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The effective field theory approach allows a rigorous disentangling of high and low 
energy effects in the heavy quarkonium dynamics. Focusing in particular on the 
spectrum, we describe the nature of the non-perturbative effects and discuss our 
present knowledge of them. 



1 Introduction 

It is well known that the complex structure of the QCD vacuum can, in some 
regimes, be parametrized by local vacuum condensates, i.e. the expectation 
values of operators where all the perturbative vacuum fluctuations are taken 
out at the best of our ability 111. The first attempt to calculate the effect of 
non-perturbative vacuum condensates on the energy levels of heavy quarkonia 
was performed in el. The leading contribution (in a s and Aqcd , the scale of 
non-perturbative physicsl of the vacuum condensates to Coulombic n, I heavy 
quarkonium states reads □: 

where e„ is a (known) number of order I and G 2 = ((a s / ' n)F^ v (Q)F a ^ '(0)) is 
the gluon condensate. The above expression for the corrections to Coulombic 
energy levels displays two relevant characteristics: 1) it is a correction of non- 
potential type, like the Lamb shift in QED; 2) it grows like ~ n 6 , thua_bcing out 
of control for levels beyond the ground state. It was soon realized Bu that the 
strong growth in n could be corrected to some extent by considering non-local 
gluon condensates 

G 2 {x) ~ (^F- v (x)<f>(x,0)^F b ^(0)) , (2) 
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where (f>(x, 0) ad -' is the Schwinger line (in the adjoint representation) connecting 
x with 0. These were understood as due to the presence of a fluctuating gluonic 
background with a characteristic time length T g ~ Aq CD . It is apparent that 
the local condensates stem from an expansion of the non-local ones in cases in 
which the correlation length is large with respect to the other physical scales 
of the system. 

Heavy quarkonium is a non-relativistic bound system. Besides Aqcd, it is 
characterized by at least three hierarchically ordered scales, mv 2 <C mv <C m, 
where m is the mass of the heavy quark and v its velocity. Therefore, it 
is only under the condition Aqcd <C mv 2 that the use of local condensates 
and, hence, the Voloshin-Leutwyler formula (||) can be justified. In the other 
regimes where heavy quarkonium states can sit, the non-perturbative dynamics 
will be encoded into more extended objects: non-local condensates and Wilson 
loop operators. This is the case for higher quarkonium levels (n > 1). 

In the following we will consider non-perturbative effects in heavy quarko- 
nium systems in different kinematic regimes, in an effective field theory context. 
This approach has not only the considerable practical advantage of disentan- 
gling the different dynamical scales of the system, but also the conceptually 
relevant feature of disentangling perturbative effects from non-perturbative 
ones. This allows us to fully exploit the predictive power of QCD. 

2 Local and non-local condensates: Aqcd "C mv 

For the lower lying quarkonium states, it can be expected that the inverse of 
the typical size of the system is larger than Aqcd- If this condition is fulfilled, 
both scales m and mv can be integrated out perturbatively from QCD, leading 
to an effective field theory where only the degrees of freedom of order Aqcd 
or mv 2 remain dynamical. This effective field theory is known as potential 
NRQCD, pNRQCDM. 

The infrared sensitivity of the quark-antiquark static potential at three 
loops 13 signals that it may become sensitive to non-perturbative effects if the 
next relevant scale after mv is Aqcd- Indeed, in the situation mv 3> Aqcd 3> 
tow 2 , the leading non-perturbative contribution (in a s and in the multipole 
expansion) to the static potential reads Era 

2 2 roc 

Vo(r) non - port = -i^-T F — / dte~ ltCAa ^^(E a (t)<j)(t,0)lf-E b (0))(p), (3) 
N c 3 Jo 

where r is the quark-antiquark distance. This term explicitly cancels, up to 
the considered order, the dependence of the perturbative static potential on the 
infrared scale /_t. It is interesting to note that the leading contribution in the 
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Aqcd /mv 2 expansion of y non -P ort (obtained by putting the exponential equal 
to 1) cancels the order Aq CD r 2 renormalon that affects the static potential 
(the leading-order renormalon, of order Aqcd , cancels against the pole mass) . 
Therefore, also in renormalon language, the above operator is the relevant non- 
perturbative contribution to the static potential in the considered kinematic 
situation. 

If Aqcd < m» 2 the static potential is purely perturbative and its explicit 
dependence on the infrared scale /x is reabsorbed in a physical observable by 
non-potential contributions. In the specific caseL of the quarkonium energy 
levels up to order a!? ln/i, these contributions areEl 

SE ntl>j = -i^-T F J dt(n,/|re it ^- ff »)r|r 1 ,0(E a (t)^(t ) 0)^E 6 (0))(M),(4) 
bj„. — — m— — — , ti n 



4n 2 ' to 2N C r ' 

where \n, I) are the Coulomb wave functions. It is worth while to notice that, 
if Aqcd *C mv 2 , the scale mv 2 can be integrated out perturbatively from the 
above formula and the non-perturbative contributions reduce to the Voloshin- 
Leutwyler formula (0), as can easily be seen by recognizing that: 



{mC F a s ) A 3N, 



Tf ( n, I 
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The above outline leads to the following conclusion. For quarkonium of 
a typical size smaller than 1/Aqcd, the most relevant operator of the non- 
perturbative dynamics is the bilocal gluon condensate (E a (t)<^(£, 0)^ j E fc (0)}, 
which belongs to the class of non-local gluon condensates considered in the 
introduction. In the following section we will discuss our present knowledge of 
it. 



2.1 The non-local condensate (F£ v (x)4>(x, 0)^^ (0)) 

The correlator (F^Jx)(j)(x, 0) ab ^F| (0)) is perturbatively known at the next-to- 
leading order in a s cB. However, here we are interested in its noi^pitjjjjhali- 
ve behaviour. Different parametrizations have been proposed E30'olijii3. 
Because of its Lorentz structure, the correlator is in general described by two 
form factors. A convenient choice of these consists in the chromoelectric and 
chromomagnetic correlators: 

(E a (x)(f)(x, 0)^ j E b (0)> , (B a (x)(j)(x, 0)^ j B b (0)}. 
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The strength of the correlators is of the order of the gluon condensate. In the 
long range (a; 2 — > oo) they fall off exponentially (in the Euclidean space) with 
some typical correlation lengths. In the following we will concentrate on these 
correlation lengths. 

The lattice calculation E3, using cooling techniques, obtains the same cor- 
relation length (T g ) for both form factors, and this is 

T g = 0.34 ± 0.02 ± 0.03 fm (4 flavours, am = 0.01), (5) 
T g = 0.22 ± 0.01 ± 0.02 fm (quenched). (6) 

The less accurate, but traditional (quenched) lattice calculation done in 
obtains two different correlation lengths for the chromoelectric {Tg) and the 
chromomagnetic correlators {T^): 

T E ^ Tf ~ 0.1-0.2 fm (quenched). (7) 

Finally, a recent sum-rule estimation obtains Tf < Tf. The sum rule 
turns out not to be stable for the chromoelectric correlator, while for the 
chromomagnetic correlation length it gives 

T g B = 0.13t8;{J|j fm (3 flavours), (8) 
T g B = O.lllo °2 fm (quenched). (9) 

The correlation lengths Tg and Tg have a precise physical interpretation. 
Their inverses correspond to the masses of the lowest-lying vector and pseu- 
dovector static quark-gluon hybrids, respectively. This can be explicitly seen 
in the short-range limit, x — > 0, where the hybrida-^in this case also called 
gluelumps) operators can be explicitly constructed BO. The suitable effective 
field theory is pNRQCD in the static limit Gluelump operators are of the 
type Tr{07?}, where O = O a T a corresponds to a quark-antiquark state in the 
adjoint representation (the octet) and H = H a T a is a gluonic operator. By 
matching the QCD static hybrid operators into pNRQCD, we get the static 
energies (also called potentials) of the gluelumps. At leading order in the 
multipole expansion, they read 

V H {r) = V Q (r) + -L, (10) 

9 

(i? a (i)(Ht,0)^ j ff fc (0)) non -P Grt - he~ lt ' T " +.... 

Since hybrids are classified in QCD according to the representations of -Doo,h, 
while in pNRQCD, where we have integrated out the length r, their classifi- 
cation is done according to the representations of 0(3) x C, the static hybrid 
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shor±=xange spectrum is expected to be more degenerate than the-lpng-range 
one E30. The lattice measure of the hybrid potentials done in Eia confirms 
this feature. InQ it has been shown that the quantum numbers attribution of 
pNRQCD to the short-range operators, and the expected 0(3) x C symmetry 
of the effective field theory match the lattice measurements. By using only E 
and B fields and keeping only the lowest-dimensional representation we may 
identify the operator H for the short-range hybrids called £+ (and II g ) with 
r • E (and r x E) and the operator H for the short-range hybrids called E~ 
(and n u ) with r • B (and r x B). Hence, the corresponding static energies for 
small r are 

V.n. W = Vo{r) + 7^' ^.n>) = V {r) + -L. 

a a 

The lattice measure of@ shows that, in the short range, V + > (r) > V s - n (r) . 

This supports the sum-rule predictionliJ that the pseudovector hybrid lies lower 
than the vector one, i.e. T® < Tg . 

3 Wilson loop operators: Aqcd ~ M 

For higher quarkonium levels, Aqcd is expected to be comparable with mv. We 
cannot match into pNRQCD perturbatively, since the scale associated to the 
quarkonium size r is already non-perturbative. The relevant non-perturbative 
dynamics is therefore contained in more extended objects than (local or non- 
local) gluon condensates: Wilson loops and field insertions on these. 

In particular, disregarding effects due to scales lower than Aqcd, the static 
potential is given byl23 



V (r) = lim ^ln(Wn), (11) 

where Wn is the static Wilson loop of size r x T and ( ) means an average 
over the gauge fields. Lattice studies tell us that at distances r ~ 1/Aqcd the 
potential is no longer Coulombic but rises linearly (Vo(r) ~ err). Higher^order 
corrections in the 1/m expansion have been calculated over the years EU and 
are given by field strength insertions on the Wilson loop. For instance the 
next-to-leading potential in the 1/m expansion isr^ 



fT/2 


,T/2 


/ dt 


/ dt'\t-t'\ 


J-T/2 J 


-T/2 



;(E(t).E(f))) n (12) 



•((E(i))) D • «E(0»n 



(13) 
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where (( ))□ means a normalized gauge average in the presence of the static 
Wilson loop. 

Wilson-loop operators of the above type may be interpreted as a piperpo- 
sition of states, describing gluonic excitations between static sourcesE3. They 
have been so far evaluated only_inside QCEL vacuum models or by lattice sim- 
ulations (for some reviews, seeEij and alsoEfl). 

3.1 The stochastic expansion 

We can now ask if some relation can be established between the Wilson loop 
(and field strength insertions on it) and the non-local gluon condensates dis- 
cussed above. In fact we can express the Wilson loop as a formal expansion 
in terms of ghionic correlation functions by means of the so-called stochastic 
expansion EIO 

ln(^ D )=V^ / d5„ 1 ^(« 1 )-.-dS^(u n )^(0,«i) 

xf 1 " 1 (uiMm, o) • • • 0(0, u n )F^ u " K)0(«„, o)) cum , (u) 

where S(d) denotes a surface whose contour the rectangular Wilson loop. The 
cumulants ( } cum are defined as 

(0(O,ui)F( Ul )0( Ml ,O)) cum = (0(O,ui)F( Ul )0( Ul ,O)) = 0, 

(0(0, Ui)F(ui)(f)(ui, U 2 )F(U 2 )(I>(U2, 0))cum = 

(0(0, u 1 )F{u x )<f>{u 1 , u 2 )F{u 2 )<t>{u2, 0)) 
-(0(0, u 1 )F(u 1 )i>{u 1 ,0)) (0(0, u 2 )F(u 2 )<f>(u 2 , 0)) . . . = 
( J F( Ml )0( Wl , U2 ) adj F( U2 )), 

It is important to realize that the expansion of Eq. ([l4|) is substantially differ- 
ent with respect to the expansions in 1/m and r, which led to the construction 
of the low-energy effective field theories discussed above. Those expansions 
were justified by the dynamics of the system under study, and each term of 
it is, indeed, suppressed by powers of the highest dynamical scale left divided 
by the scale (which is larger) that has been integrated out. Instead, from a 
power counting point of view, each term of the expansion ( |l4| ) is of the same 
size (for instance each term of the series is in general expected to contribute 
to the string tension a). As soon as we truncate the series, e.g. up to the 
bilocal cumulant (the first non- vanishing one), we introduce an uncontrolled 
approximation and define a model. This model is known as the model of the 
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stochastic vacuuniLS. Indeed, this model turns out to be quite successful in the 
study of processes that involve quark systems that may be described by almost 
static Wilson loops (for some reviews, seeo). Let us only mention that the 
model predicts, in agreement with the lattice data, a long-range linear static 
potential with slope a ~ T® 2 Gi. It would be highly desirable to have a field 
theoretical justification of the expansion (|l4|); however, such a justification is 
missing up to now (for a recent investigation on higher cumulants, seec3). 

4 Conclusions 

In the previous sections we have shown how the non-perturbative QCD vacuum 
enters the dynamics of heavy quarkonium. As much as the non-perturbative 
scale Aqcd is bigger than the dynamical scales of the non-relativistic system 
(to, mv and mv 2 ), as extended the relevant non-perturbative operators are. 
In the situation Aqcd <C mv 2 these operators reduce to local condensates, 
i.e. some numbers. This is the most favourable situation for a theoretical 
investigation. The bottomonium and charmonium gmund states have been 
investigated in this framework (for a recent review, seeEH). In the limiting case 
of ti threshold production the non-perturbative corrections can be -neglected 
and the investigation is completely accessible to perturbative QCDcS. In the 
situation mv 2 < Aqcd mv the non-perturbative physics is encoded into 
non-local condensates. The dominant one is the bilocal gluon condensate. 
As we have discussed above it essentially depends on its strength, i.e. the 
gluon condensate, and on some correlation lengths, which can be related to the 
masses of the lowest-lying (static) quark-gluon hybrid resonances. Although 
this situation seems quite interesting, it has been poorly investigated in-heavy 
quarkonium phenomenology, and mainly in the framework of models Eru. The 
main reason is, in our opinion, that, while the quarkonium ground state seems 
to be accessible by a purely perturbative treatment (plus local condensates) 
and the potential describing higher excited quarkonium states seems entirely 
dominated by non-perturbative effects, it is, a priori, not clear to which states 
the situation mv 2 < Aqcd <C mv is applicable. We will come back to this 
point at the end of this section. Finally in the situation Aqcd ~ mv the 
relevant non-perturbative objects are Wilson-loop operators. A lattice study 
of them and the subsequent quarkonium spectroscopy have been done incj. 

The outlined study is rigorous and allows a systematic disentanglement 
of the high from the low energy scales of the heavy quarkonium system under 
study. In this specific sense also perturbative and non-perturbative effects turn 
out to be disentangled. 

We conclude by mentioning a somehow delicate point, connected with the 
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scales of heavy quarkonium systems: the difficulty to state a priori in which 
kinematic situation a particular quarkonium state is. This is mainly due to the 
fact that the scales mv and mv 2 are not so well defined, nor so widely separated 
that different kinematic situations cannot overlap in a preliminary analysis Ifpr 
a related discussion on the energy scales also in quarkonium physics, seeEJ). 
Therefore, there are situations in which scales can only be fixed a posteriori, 
i.e. by assuming a particular situation and by checking that the final result is 
consistent with it and with the experimental data. 
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